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A general approach to analyze the electrodynamics of nuclear matter in bulk is presented 
using the rclativistic Thomas- Fermi equation generalizing to the case of ~ (jripianck/^T-n)'^ 
nucleons of mass m„ the approach well tested in very heavy nuclei (Z ~ 10^). Particular 
' attention is given to implement the condition of charge neutrality globally on the entire 
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configuration, versus the one usually adopted on a microscopic scale. As the limit N 



, (wpianck/'Tin)'^ is approachcd the penetration of electrons inside the core increases and a 

^ '• relatively small tail of electrons persists leading to a significant electron density outside 
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the core. Within a region of 10^ electron Compton wavelength near the core surface electric 
fields close to the critical value for pair creation by vacuum polarization effect develop. These 
results can have important consequences on the understanding of physical process in neutron 
stars structures as well as on the initial conditions leading to the process of gravitational 
collapse to a black hole. 



Qh| It is well know that the Thomas-Fermi equation is the exact theory for atoms, molecules and 

solids as Z ^ oo 1] . We show in this letter that the relativistic Thomas- Fermi theory developed for 



o 

^ I the study of atoms for heavy nuclei with Z ~ 10^ gives important basic new information 

on the study of nuclear matter in bulk in the limit of ~ ("^Planck /""^n)^ nucleons of mass 
^ I run and on its electrodynamic properties. The analysis of nuclear matter bulk in neutron stars 



composed of degenerate gas of neutrons, protons and electrons, has traditionally been approached 
by implementing microscopically the charge neutrality condition by requiring the electron density 
rieix) to coincide with the proton density np{x), 

ne{x) = np{x). (1) 

It is clear however that especially when conditions close to the gravitational collapse occur, there 
is an ultra-relativistic component of degenerate electrons whose confinement requires the existence 
of very strong electromagnetic fields, in order to guarantee the overall charge neutrality of the 
neutron star. Under these conditions equation will be necessarily violated. We are going to 
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show in this letter that they will develop electric fields close to the critical value Ec introduced by 
Sauter Q|, Heisenberg and Euler Q], and by Schwinger Q| 
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Special attention for the existence of critical electric fields and the possible condition for electron- 
positron (e'*"e~) pair creation out of the vacuum in the case of heavy bare nuclei, with the atomic 
number Z > 173, has been given by Pomeranchuk and Smorodinsky [21, Gershtein and Zel'dovich 

in ^ n 

J], Popov Popov and Zel'dovich 6], Greenberg and Greiner 'IQ\, Muller, Peitz, Rafelski and 
Greiner They analyzed the specific pair creation process of an electron-positron pair around 
both a point-like and extended bare nucleus by direct integration of Dirac equation. These con- 
siderations have been extrapolated to much heavier nuclei Z ^ 1600, implying the creation of a 
large number of e'^e'^ pairs, by using a statistical approach based on the relativistic Thomas- Fermi 
equation by Muller and Rafelski 8], Migdal, Voskresenskii and Popov 9]. Using substantially the 
same statistical approach based on the relativistic Thomas- Fermi equation, Ferreirinho et al. lll |. 
Ruffini and Stella |l^ have analyzed the electron densities around an extended nucleus in a neutral 
atom all the way up to Z ~ 6000. They have shown the effect of penetration of the electron orbitals 
well inside the nucleus, leading to a screening of the nuclei positive charge and to the concept of an 
"effective" nuclear charge distribution. All the above works assumed for the radius of the extended 
nude, the .....PMC.: f„..u>a.Q, 

Rc w roA^/^ ro = 1.2 • 10"^^cm, (3) 

where the mass number A = Nn + -/Vp, Nn and Np are the neutron and proton numbers. The 
approximate relation between A and the atomic number Z = Np, 

Z^^, (4) 
was adopted in Refs. 0,0], or the empirical formulae 

^-4 + 25oI^73l"'> 



was adopted in Refs. 
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121. 



The aim of this letter is to outline an alternative approach of the description of nuclear matter 
in bulk: it generalizes, to the case of ~ (jn-pianck/^n)^ nucleons, the above treatments, already 
developed and tested for the study of heavy nuclei. This more general approach differs in many 
aspects from the ones in the current literature and recovers, in the limiting case of A smaller than 
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10^, the above treatments. We shall look for a solution implementing the condition of overall 
charge neutrality of the star as given by 

Ne = Np, (6) 

which significantly modifies Eq. since now Ne{Np) is the total number of electrons (protons) 
of the equilibrium configuration. Here we present only a simplified prototype of this approach. 
We outline the essential relative role of the four fundamental interactions present in the neutron 
star physics: the gravitational, weak, strong and electromagnetic interactions. In addition, we 
also implement the fundamental role of Fermi-Dirac statistics and the phase space blocking due 
to the Pauli principle in the degenerate configuration. The new results essentially depend from 
the coordinated action of the five above theoretical components and cannot be obtained if any one 
of them is neglected. Let us first recall the role of gravity. In the case of neutron stars, unlike 
in the case of nuclei where its effects can be neglected, gravitation has the fundamental role of 
defining the basic parameters of the equilibrium configuration. As pointed out by Gamow [l^ . 
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in general relativity, configurations of 



at a Newtonian level and by Oppenheimer and Volkoff 
equilibrium exist at approximately one solar mass and at an average density around the nuclear 
density. This result is obtainable considering only the gravitational interaction of a system of 
Fermi degenerate self-gravitating neutrons, neglecting all other particles and interactions. It can 
be formulated within a Thomas-Fermi self-gravitating model (see e.g. In the present case 

of our simplified prototype model directed at evidencing new electrodynamic properties, the role 
of gravity is simply taken into account by considering, in line with the generalization of the above 
results, a mass-radius relation for the baryonic core 

rNS ^ ^ ^ rnpi^ncl. ^ 
m,rC ran 

This formula generalizes the one given by Eq. ^ extending its validity to ~ (w-pianck/'^T'n)^) 



leading to a baryonic core radius Rc ~ 10km. We also reca^ 
nuclear matter in bulk in neutron stars ( see e.g. Bethe et al. 
at mass densities larger than the "melting" density of 



I that a more detailed analysis of 



3 and Cameron Q 



) shows that 



Pc = 4.34 • lO^'^y/cm^ (8) 

all nuclei disappear. In the description of nuclear matter in bulk we have to consider then the 
three Fermi degenerate gas of neutrons, protons and electrons. In turn this naturally leads to 
consider the role of strong and weak interactions among the nucleons. In the nucleus, the role of 
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the strong and weak interaction, with a short range of one Fermi, is to bind the nucleons, with a 
binding energy of 8 MeV, in order to balance the Coulomb repulsion of the protons. In the neutron 
star case we have seen that the neutrons confinement is due to gravity. We still assume that an 
essential role of the strong interactions is to balance the effective Coulomb repulsion due to the 
protons, partly screened by the electrons distribution inside the neutron star core. We shall verify, 
for self-consistency, the validity of this assumption on the final equilibrium solution we are going 
to obtain. We now turn to the essential weak interaction role in establishing the relative balance 
between neutrons, protons and electrons via the direct and inverse /3-decay 

p + e — > n + Ve, (9) 
n — > p + e + Uf,. (10) 

Since neutrinos escape from the star and the Fermi energy of the electrons is null, as we will show 
below, the only non-vanishing terms in the equilibrium condition given by the weak interactions 
are: 

[(P„^c)2 + Mlcr' - M^c' = [{P^cf + Mlc'fl' - M,^ + \e\V^^^„ (11) 

where and Pp are respectively, the neutron and proton Fermi momenta, and V^o^l is the 
Coulomb potential of protons. At this point, having fixed all these physical constraints, the main 
task is to find the electrons distributions fulfilling in addition to the Dirac-Fermi statistics also 
the Maxwell equations for the electrostatic. The condition of equilibrium of the Fermi degenerate 
electrons implies the null value of the Fermi energy: 

[(Pf c)2 + m^c^] 1/2 _ ^ ^Vconiir) = 0, (12) 
where Pf is the electron Fermi momentum and V^,niii (^) the Coulomb potential. In line with the 
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procedure already followed for the heavy atoms 
Fermi Equation: 

Tn— = -47ra < 9{x - x^) - 



12l | we here adopt the relativistic Thomas- 



2 n3/2^ 



(13) 



where a = /{he), 9{x — Xc) represents the normalized proton density distribution, the variables 
X and X ^'^^ related to the radial coordinate and the electron Coulomb potential V^oui by 

' ' ; eVcoui{r) = , (14) 



Rc V 47r 
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and the constants Xc{r = Rc) and f3 are respectively 



4,7 • ^^^{w;} ■ ('^) 

The solution has the boundary conditions 

X(0) = 0; x(oo) = 0, (16) 

with the continuity of the function x s-i^d its first derivative x' at the boundary of the core Rc- 
The crucial point is the determination of the eigenvalue of the first derivative at the center 

X'(0) = const, (17) 

which has to be determined by fulfilling the above boundary conditions (jlGj) and constraints given 
by Eq. (|11|) and Eq. 0. The difficulty of the integration of the Thomas- Fermi Equations is 
certainly one of the most celebrated chapters in theoretical physics and mathematical physics, 
still challenging a proof of the existence and uniqueness of the solution and strenuously avoiding 
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the occurrence of exact analytic solutions. We recall after the original papers of Thomas 

{2^, the works of Scorza Dragoni j^l; Sommerfeld 25 1, Miranda 2^ all the way to the many 
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Fermi |23l |. the works of Scorza Dragoni [2J|, Sommerfeld 

hundredth papers reviewed in the classical articles of Lieb and Simon []J, Lieb |2|3] and Spruch 
The situation here is more difficult since we are working on the special relativistic generalization 
of the Thomas-Fermi Equation. Also in this case, therefore, we have to proceed by numerical 
integration. The difficulty of this numerical task is further enhanced by a consistency check in 
order to fulfill all different constraints. It is so that we start the computations by assuming a total 
number of protons and a value of the core radius Rc- We integrate the Thomas- Fermi Equation 
and we determine the number of neutrons from the Eq. Hll|) . We iterate the procedure until a 
value of A is reached consistent with our choice of the core radius. The paramount difficulty of the 
problem is the numerical determination of the eigenvalue in Eq. (|17|) which already for A ~ 10^ had 



presented remarkable numerical difficulties 1]|. In the present context we have been faced for a 
few months by an apparently unsurmountable numerical task: the determination of the eigenvalue 
seemed to necessitate a significant number of decimals in the first derivative (|17|) comparable to 
the number of the electrons in the problem! We shall discuss elsewhere the way we overcame the 
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difficulty by splitting the problem on the ground of the physical interpretation of the solution 
The solution is given in Fig. and Fig. ((21). 

A relevant quantity for exploring the physical significance of the solution is given by the number 
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FIG. 1: The solution x of the relativistic Thomas- Fermi Equation for A — 10^^ and core radius Rc — 10km, 
is plotted as a function of radial coordinate. The left red line corresponds to the internal solution and it is 
plotted as a function of radial coordinate in unit of Rc in logarithmic scale. The right blue line corresponds 
to the solution external to the core and it is plotted as function of the distance Ar from the surface in the 
logarithmic scale in centimeter. 

of electrons within a given radius r: 

Ne{r) = r ATr{r'fne{r')dr' . (18) 
Jo 

This allows to determine, for selected values of the A parameter, the distribution of the electrons 
within and outside the core and follow the progressive penetration of the electrons in the core at 
increasing values of A [ see Fig. (j^J]. We can then evaluate, generalizing the results in 
the net charge inside the core 



nl, 
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N^et = Np - Ne{Rc) < Np, (19) 

and consequently determine of the electric field at the core surface, as well as within and outside 
the core [see Fig. @] and evaluate as well the Fermi degenerate electron distribution outside the 
core [see Fig. ((SJ]. It is interesting to explore the solution of the problem under the same conditions 
and constraints imposed by the fundamental interactions and the quantum statistics and imposing 
instead of Eq. ^ the corresponding Eq. Q. Indeed a solution exist and is much simpler 

Unix) = np{x) = ne{x) =0, X = 0. (20) 
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FIG. 2: The same as Fig. enlargement around the core radius Rc showing exphcitly the continuity of 
function x its derivative x' from the internal to the external solution. 

Before concluding as we announce we like to check on the theoretical consistency of the solution. 
We obtain an overall neutral configuration for the nuclear matter in bulk, with a positively charged 
baryonic core with 

A.„.=„..2(™)^(^)Vi)\ (21) 

and an electric field on the baryonic core surface (see Fig. Q) ) 

-| = 0.92. (22) 

The corresponding Coulomb repulsive energy per nucleon is given by 

U^^ = ^ {^y^^' « 1-78 • lO-^(MeV), (23) 

well below the nucleon binding energy per nucleon. It is also important to verify that this charge 
core is gravitationally stable. We have in fact 

Q_ = 1.56 . 10-4. (24) 

VGM \m^J ^ ^ 

The electric field of the baryonic core is screened to infinity by an electron distribution given 
in Fig. ((21). As usual any new solution of Thomas-Fermi systems has relevance and finds its 
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FIG. 3: The electron number (|18|l in the unit of the total proton number Np, for selected values of A, is 
given as function of radial distance in the unit of the core radius i?c, again in logarithmic scale. It is clear 
how by increasing the value of A the penetration of electrons inside the core increases. The detail shown in 
Fig. Q and Fig. |(SJ demonstrates how for N ~ (TOpianck/"^n)^ a relatively small tail of electron outside the 
core exists and generates on the baryonic core surface an electric field close to the critical value given in . 
A significant electron density outside the core is found. 

justification in the theoretical physics and mathematical physics domain. We expect that as in 
the other solutions previously obtained in the literature of the relativistic Thomas-Fermi equations 
also this one we present in this letter will find important applications in physics and astrophysics. 
There are a variety of new effects that such a generalized approach naturally leads to: (1) the mass- 
radius relation of neutron star may be affected; (2) the electro dynamic aspects of neutron stars and 
pulsars will be different; (3) we expect also important consequence in the initial conditions in the 
physics of gravitational collapse of the baryonic core as soon as the critical mass for gravitational 
collapse to a black hole is reached. The consequent collapse to a black hole will have very different 
energetics properties. 
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